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CHAPTER 1 


INTRODUCTION 


Turbomachine blade vibration has been an important tedinical field for many years 
Blades m a turbomadime represent the most cntical componaits m the design of rotating 
machmery With mcreasing seventy m operating aivironmait of fixture turbomadhines, the need 
for more effiaent blade design is essaidal 

A typical blade is pretwisted along its laigth, has curvature (or camber) in the diordwise 
direction and vanes m thickness along the chord. TTie cross-sections are generally aerofoil m 
shape Both the thickness and the curvature vanes along the length, but decrease from the root to 
the end The blade width also dianges along the laigth. 


Vibration induced fatigue failure of blades is a problem of major concern in turbomachine 
bladmg The bdiaviour of most blades is strongly influenced by the mean and alternating stresses 
expenenced during the operation, by the applied loading history, by the elastic and fatigue 
properties of the blade and by the operating environmait. The cntical aspect of the problem 
constitutes resonant vibration at an mt^al order speed when the principal conponent of the 
vibratmg modes in the spinning blade matdies both m time and space a correspondmg pattern in 
the gas/steam stream , caused by any distortion m the flow field. This may give nse to large 
dynamic stresses leadmg to fatigue failure. The practice of tunmg the blade away from natural 
frequencies may not always be possible in an engine of advanced design with as many as thousand 
blades of different diaractenstics As compared to typical steam turbine power generating 
systems, with heavy rotors and low operating ranges of speeds, the nature of the problem become 
more apparent for hi^ performance gas turbmes operating in superscmic regimes containing 
several cntical speeds The blades get excessively stressed not only during operation of the rotor at 
constant speeds, but also dunng their transit past several criticals during acceleration and 
deceleration operations of the turbomachine In this study, the problem of the drtermmation of 
transiait resonant stresses, during variable speed operations like start-up and shut-down of the 
madiine, is addressed 
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available unsteady force components and checked his results experimentally. Hoyniak and Fleeter 
[1981] made an energy balance between unsteady aero-dynamic work and energy dissipation due 
to aerodynamic damping to predict blade resonant vibrations by employing experimental unsteady 
aerodynamic gust data for flat chamber cascade airfoils Jadwani and Rao [1982], analyzed the 
force vibration response using the equations of Lagrange They designed an expenmaital rig 
[1982] simulating nozzle by permanent magnets and measured the tip deflection of the rotating 
blades under nozzle excitation Hinyuki Kojima [1985] considered the case of forced vibration of 
a beam with a mass subjected to altematmg electro-magnetic force and solved the governing 
partial differential equatons by harmonic balance method Leung [1983] used the dynamic 
stiffeess method relating the amplitude of applied forces and response of a harmonically vibrating 
contmuum dampmg was mcorporated along with load distribution and equivalence of the solution 
to Duhamel's mtegral was illustrated. Comparison was made with modal analysis m^od. 

Mac Bam [1976] used finite differaice method to determine the transient response of a 
cantilever plate subjected to normal impact by a ballistic paidulum Sisto and Chang [1981] 
considered the gyroscopic forces on the dynamic bdiavior of rotatmg blades, Coriolis forces were 
mcluded m the analysis of Sisto, Chang and Sutcu [1982].They also mvestigated the instabilities 
due to harmonic vanation with time of the processional rate due to whirling and other causes 
Irretier [1984], usmg beam theory and modal approximation technique, presetted a numerical 
modal for an untwisted blade to obtain transient vibration response due to partial admission durmg 
a change of rotational speed Banneijee and Kennedy [1985] derived the theoretical expressions for 
the displacement response of an axially loaded beam subjected to concentrated or distributed loads 
having stationary and ergodic properties 

While considerable amount of activity is witnessed m obtaining the forced vibration 
response of blades under constant operational speeds of the rotor, transiait resonant stress analysis 
durmg variable speed operations has not beai able to receive any sigmficant attaition Sudi an 
analysis is of considerable inportance for aircraft turbine bladmg, which experience frequent 
dianges m rotor speed, durmg various manoevours of the aircraft 

1.2 Present Work 

In the present study an attempt has been made to develop a procedure for estimation of 
turbme blade stresses during typical variable speed turbme operations like start-up and shut down. 
Initially, the r^ponse of a single degree of freedom mass-spring-dashpot model, to an external 
force, with a time depaideit exdtation fi-equaicy, is obtained. For ease of analysis, the exdtation 
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frequency is taken to vary uniformly with time The response is obtained through the use of 
Duhamel's Integral The response is obtained for different values of damping ratios for vanous 
frequency sweep rates, in non-dimensional form The analysis is then extended to the case of an 
actual turbine blade Using beam theory, n^ecting hi^er order eiffects such as shear deformation, 
rotary mertia, Conolis forces etc , the lonetic energy and the Reissnef s functional for a typical 
tapered, twisted, asymmetnc aerofoil cross-section blade mounted at a stagger angle on a disc 
under rotational acceleraton have been denved. Available penodic excitation data , simulating 
nozzle exatation, is used m Founer form m the traction forces term of the Reissnef s functional 
Ritz process has been employed to obtam the equations of forced motion. The eigenvalue problem 
IS solved to obtam the free vibration characteristics Modal Analysis tedinique is then used to 
decouple the equations of motion The nondimensional response obtained for the single d^ee of 
freedom system, is processed with the modal equations of the blade to obtam the transient resonant 
strresses inflicted on the blade dunng the vanable speed operation of the turbomachine Results 
are obtamed for typical acceleration rates and compared with the correspondmg resonant stresses 
for constant speed operations 



CHAPTER 2 


RESPONSE OF A SINGLE DEGREE FREEDOM SYSTEM 
WITH VARIABLE FREQUENCY EXCITATION FORCE 


The response of a single d^ee freedom, mass-spnng-dashpot system, to an excitation 
force with frequaicy varymg with Pme is analysed, m this chapter. The vanation in the excitation 
frequaicy is taken to be constant with time, for ease of analysis. Such an analysis is essential, 
pnor to analysmg the turbme blade - firstly to understand the phaiomenon of a momentanly 
expenenced resonance and secondly, to overcome the computational problems, descnbed m the 
next chapter, posed by the turbme blade calculations 


2.1 Equation of motion and response 


Consider a smgle d^ee freedom system as shown m Fig 2 1 with an exatation force, 
jp(t), such that at any instant of time, the mstantaneous excitation frequency, (o , can be 
expressed as 


<o = 6)o +ai 

The mstantaneous force is 


¥{{) = jFq sm[ m 0 ^ + 2 


(2 1 ) 


(2 2 ) 


<y 5 is the initial excitation frequency, a is the sweep rate or acceleration. The magnitude 
Fg, of the excitation force is assumed to be constant. (A cosine excitation c:an be readily 
accommodated m the above expressions by includmg appropnate phase angle within the sine 
excitation force term.) 


The equation of motion can be written as 
mx + cx + kx = Fo sin|^6Jo/ + -ia/^j . 


(2 3) 


The above equation of motion does not have an exact analytical solution Referoice can be made 
to Baker [1939], for a graphical method of obtaining the response A numerical procedure is 
adopted in the present study. 
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Fig. 2.1 


Single degree of freedom system with variable frequency excitation. 
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Defining the following non-dimensional parameters 

T = pt = non - dimensional time 

Tq = = initial frequency ratio 

P 


r^- — = instantaneous frequency ratio 
P 

= nondimensional sweep rate 
P 


x(t) 

H = = dynamic magmfication factor (or displacenent ratio) 

p 

^ = static deflection 
mp 


equation (2 3) can be written as 


d^H ^ dH , 1 2 , 

-^ + 2f— +« = sm(r,r+-r.i- ) 

For zero initial conditions, the response of the system, usmg Duhamel's Integral, is 


/f(r) = |e 
0 




sinVr^^(r-r')sin[^ror'+ir„i:'^J dt' 


(2 4) 


(2 5) 


(2 6 ) 


For 1 the above becomes 

= sin(r- r')sin^ror' <fr' (2.7) 

2.2 Numerical results 

Ten-pomt Gaussian-quatrature is employed to evaluate the mtegral in equation (2 7) and 
the nondimensional response (displacemait ratio) is obtained as a function of nondimensional time. 
Figures 2.2 - 2 5 show the nondimensional response H as z function of the nondimensional 
mstantaneous excitation frequaicy, r, (= + r„r) , for various values of sweep rate, ,and the 

dampmg ratio, ^ . 

As is known, for a case whai exataticm frequency is constant (i.e. a = 0), resonance 
occurs vshai co = p hi the absaice of damping, infinite amplitude of vibration can be 
theoretically obtained whai {ojp = 10 However it takes time for these amplitudes to build iqi As 
observed from Fig.2.2(a), the amplitude of vibration builds up as the mstantaneous frequency 
Kf approaches 1.0 . At = 1.0 infinitely large amplitude would have built up, had the system 
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been allowed to 'dwell' at the point The amplitude of vibration starts building ip at this point, but 
once past r, =10, the post resonance effect comes into play Under these influences, a large but 
finite amplitude of vibration is obtained, not at r, = 1 0 , but at the instance 'shifted' from = 1.0 
(e g at Tf -\ 06 for = 0 75 x 10'^ ) 

Figure 2 2 depicts the influence of dartipmg ratio on the peak response With a finite 
amount of dampmg existing m the system a portion of the eneigy supplied is dissipated, leaving 
less aiergy for building iqi of vibrational amplitudes and as expected 'resonant' peaks become 
lower as the dampmg mcreases 

The influaice of the sweep rate on the peak amplitude is shown in Figures 2.3. These 
figures have bean obtained for a constant value of the damping ratio (^ = 0 0025) It can be 
observed that for larger values of ,the peak amplitude get reduced and are shifted further from 
r, = 1 0 For larger , the transit bemg faster, the energy supplied is utilized more m 
acceleratmg the system than m building up vibrational amplitudes Figs 2.4 and 2 5 show similar 
results for n^ative rates of acceleration 

The peak responses for various positive and native sweqi rates, for various dampmg 
values, mcludmg those from Figs. 2.2 - 2.5 have been compiled in Table 2 1 
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Response under acceleration- Influence of the rate of acceleration 
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2.3 Remarks 

The response data generated in this chapter forms a cruaal computational con^onent for 
obtaining the blade response, described m the next chapter The nondimensional response data of 
this chapter is employed as a built-in data surface, m the blade program, whereby the blade modal 
responses are directly extrapolated from this surface 



CHAPTERS 


TURBINE BLADE RESPONSE UNDER ACCELERATION 

Tbs chapter deals with the forced vibration of a turbine blade, when the turbomachine 
speed IS mcreased or decreased , as in the step-iqi and stiqj-down operations The present analysis 
IS idealized to a case \\4iere the rate of diange m the rotational frequaicy of the turbme is constant 
Reissnefs functional, m conjunction with Ritz process has been enqiloyed to denve the equations 
of forced motion The eigaivalue problem is solved to obtain the modal properties and Modal 
analysis is subsequently earned out to decoiqile the equations of motion. The results of the 
previous chapter are then used for a computationally efficient calculation of the modal peak 
responses The stress response of the blade is conqiuted using Strength of Matenal formulae. 

3.1 Equations of motion 

Consider a typical tapered, twisted, asymmetric aerofoil cross-section blade mounted at a 
stagger angle on a disc rotating with constant acceleration shown m Figs 3 1 and 3.2. The 
centroid of the cross-section is located at G, wble O is the carter of flexure. Noting that the 
gareral displacemait of point P consists of translations about the xx and yy axes and rotation 
about the centre of flexure, O, the displacemaits field can be writtai as 

Uj, =x -yd = x^ -y 9 
1 

Uy =y + xe = yi+x 9 

, . ' (3-1) 

M, =-x u'^ -y Uy +w 
1 1 

= -x (x[-y 9')-y 9') + u 

1 .1 .1 ' 

The displacemait vector U of the particle P can be writtai as 

U = PP' = + Uyly^ + u,e, (3.2) 

\^dlere ,e^ denote the rotating umt vectors m Xj ,jPj and z directions respectively The 

unit vectors and Cy^ are related to the rotating unit vectors and in the ^ and t] 
directions as 

e =e^cJOS(j)-e^sm.<j) 
e =e^ sin^-f-c,^ cosjl 


(3.3) 



h 



Fig. 3.1 


Blade mounted on a rotating disc. 
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Fig. 3.3 Rotating vectors. 
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Equation (3 2) can be wnttai as 

= +w,e, (3 4) 

where the components of displacements and in the ^ and ;; directions are 

=(Xi ->>, ^)cos^J+(>'j +X] ^)sin^ 

(3 5) 

= -(Xj - 0 ) sm +(jj + Xj cos<^ 


The position vector of the displaced particle P' , from Figs 3 1 and 3 2 is 
F = (Re,)+(^e^ +Tje^ +^^)+(w^e^ +u^e^ +uj^) (3 6) 


The above equation can be wntten as 
F = (^+tt^)e^ +(7j+u^)e^ +(R + z + u^)e^ 

Therefore the velocity of the particle is 
P = +r^t^ +r,e, +r,t, 

Notmg from Fig 3.3 that 

e, =0)6^ =(ft)o +cct)e^ 

=-coe^ =-io)(, +at)e, 
e^=0 

equation (3 8) becomes 
P = r^e^ +ir^ +cor,)e^ +(r, -(or^)e^ 


(3 7) 


(3.8) 


(3.9) 


(3 10) 


Hence, 

= r-f 

. . . (3-11) 

= {r'l +r^ +r/)+2(<ao +at){r^i-,-r^r,)M.co^ +at) (r, +r, ) 

With the help of equations (3 7), (3.5), (3.1), the first groiqi of terms in equation (3.11) 
simplifies to 
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+y])+{A+yl)e' +2(x, j, -j, Xj)^+«" 

, , , , (3 12) 

+ Xj Xj 4'3 '’i y\ + 2Xj "“ 2 w(Xj Xj +3^j ) 

Ignoring axial inertia terms m the above 

+rl +r/ = xl + jf +{xl + yl)e'^ (3 13) 


The second term in equation (3 11) rgiresaits Conolis componaits and are neglected 
here The third group of terms m equation (3 1 1) rqiresents the kmetic aiergy due to rotation 
whidi is 


(61 0 +cctf {r^ +r/) 
= (6io +atf 


77+ (jpj + Xj ff)oos^ - (xj - jj 6 ) sin (j) 


+UR + z)u- Xi (x ' - Ti <?') - Ti iy[ + ^') 


(3 14) 


Accounting for c«itnfugal components due to displacements of the particles on the Ime of 
centroid only, the above equation simplifies to 

(61 0 +atf (r^ +r,^) 

r 1 15) 

= (6i(, +at)^ j2(/? + z)w + yf cos^ ^ + Xj sin^ (j>-x.^y^ sin2^l 


Uhaccounting for ngjd body motion and neglecting terms, with the help of equations 
(3 13) and (3 15) the kinetic energy of the blade can be written as 


T dAcb 

1 fi 1 . 

= f -pA{xi +yy) + -l^e^ +(6io +atf iR + z)upA 
o>-^ 2 

+ -(61 + at)^ pA{yi cos^ (j) + xf sin^ (f) - Xjjj sin 2^) 


dz 


The total axial displacement, w at z due to bending is given by 


(3.16) 



(3.17) 


Hence, with the help of equations (3.17) and (3.1) 
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T=p\ +[y-r,e)"]dz + ^l 




j |(x' + r'^e + r^O'^ + (j' - KO - r^ 
0 ^ 

/ 

+ sin^ (^j a{x + dz + cos^ A(y - r^6>y dz - sin 


0'Y^(R+z)dz 

t 

2f)| a{x + KyO^y - r^6)dz 

0 


(3 18) 

Using notation 

co = cOo+at (3 19) 


the kinetic aiergy can be wntten as 


+ i; 06 )^|^-}{(x' +r;0 + r0'Y +(j' -r'^-r ^')'|(;? + z)t/z 
+ sin^ dz + cos^ <j)^^A{^ - dz -sin 2<j>]A[x + r^0){y - r ^)rfzj 

(3 20) 


Reissner's Functional, which involves both the potential aiergy and complemoitary oiergy 
considerations and is known to yield simultaneously good stress and displacement fields, is 
employed along with the kmetic aiergy expression above to obtain the equations of motion for the 
blade The Reissnef s functional is givai by, [Dym and Shames (1973)], 

" III )] ^ “ Iff ^ ~ JJ ^ 

V vs 


where 


T^Sfj IS the potential energy, is the conqilementary energy, are the body 

and traction force terms respectively, w, is the virtual displacement and v and s are the body 
volume and surface respectively. 


The potential aieigy term, m the Reissner's functional is worked out to be 


v-tNlitAL 
I. I. T.. 


UiDItAit ' 
KANPUR 




i ; V ' 

«»* » m 
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j\\ryeydv = jj[T,,e,, +2 {t^£^ +T^£,J]dA± 

V- 0 A 

I 

= li-M^.x"-My + T,0")ct 
0 


(3.22) 


The complementary energy term is 






ti 


+■ 


/ 


2E 2KG 
•2 






2E{I 


yiy\ ■^*1^'! 


C.) 


2C 


(3 23) 


The body forces are assumed to be n^igible i e 


|JjA“A = 0- 


(3.24) 


Let the exatation forces on the blade be defined as components (z,t) , F^ (z, t) in the x 

and j directions and twisting moment M(z,t). These forces are periodic with the nozzle passing 
ffequaicy ff be the number of nozzles and (O the rotational frequency, the nozzle passing 
frequency is defined as 

v = n^(a = n^{a>Q+at) (3 25) 

The force componaits can be expressed m the form of Fourier senes as 


F^{x,t) = a^^{z) + Y,a„A^)cos. + +2^6„,(z)sin 

ml. “* 

Py iy, t) = (z) + 2 


, 1 2 
mn^\ CL>^t + — at 


, 1 2 
fm^\ co^t 


M(r,0 = aoM(^) + ji^(z)coa 


. 1 2 
fftfjjl (Oyt +—at 


m L 

+E*-. (z)sinj mn^ ^ 


(3 26) 
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where 


/ V' 

^Vv 

= ~ \ Fj{z,t)cx>& + dt 

y f \ 

bmxi^) = ~ J ^,(^,^)sin + — cr/^j dt 


(3.27) 


Employing equations (3 20)-(3 24) and equation (3 26) the dynamic Rassner's functional can be 


wntten as 


L = T-h 


v^pco^ -^^{{x' + r^6' + r;ef + iy' -r^0' -r^efYm,+I^)dz 

/ / ^ 

+ sin^ <j)^^ A{x + r^d)^ dz + co^ (l>^^A{y - r^6)^ dz]-%iTi2(l)^ a{^ + 


M^r, +Mlh^ Ti 

MX +My -Tg0" ^L_2L3>L_i^ tt 

- -/") 2C 

' y\y\ ^iy\ ' j 


+ [ao,(z)x + ^a„,(z)cos Tnn^[a^t + ^caA x+26„,(z)sin mnja^t + ^caA x 
m L J m L 

+ aay{z)y + Y,‘^myi^)^^ + y + 2 ] 6 „_,( 2 )sm rai.j^aiot + |at^ j y 

m L. J m L. • 

+ «0ju(^)<5' + 2!!"m^r(^)cos mn,[a}at 5>+5]i„j^(z)sin mn^{^cOf,t +^at^'^ 9)fiz 

Iff L. -J fit L. ~ 


(3.28) 


Usmg shape functions 

x=XA(0/,(^) m,=5]aW(^) 

>’ = i;s,(0/,(z) M^=2£^,(0ft,(z) (3.29) 

^=2C,(07,(z) T; =^FX)h.{z) 
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where 


J (Z) = ^ ^ ~ ^) 2 ^■*■2 1 ) 2’/+3 

' 6 3 6 

7,(Z)=z' — ^-z'"' 

/ + 1 

I + 1 

z = ^ 

/ 


which satisfy the boundary conditions of a cantilever, 
Ritz process 
dL 


dA, 


■ = 0 etc 


IS applied to the dynamic Reissner functional above to obtam the equations of motion 

[M]{9} + [C]{9} + [/q{9} = {0o + 5]{e„}cos/w«,(6)o^ + J«^' 

OT 

+5] {6m+6 )smmn^[a^t + ~at^ 


(3 30) 


(3 31) 


(3 32) 


The matnx elements have been descnbed in the Appendix 

The free vibration equation is solved computationally to obtam the natural frequenaes and 
the eigaivector matnx [U]. 

3.2 Response 

Modal analysis technique is employed to obtain forced vibration response 
Usmg 

{q}=[U]{h} (3 33) 


and premultiplymg equation (3.32) by \UY one gets 


[M]{h} + [C] W + [in W = {0OA- + S cos/rni + jot 




(3.34) 
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where 


m = [uV[M][u] 

[C] = [t/]"[C][t/] 

[X] = [uf[/q[f/] 

{Q}on =[Uf{Q}o 

{QJ, =WY{QJ 

=[f/f {2™^} 


The decoupled Ath modal equation can be readily writtai from the above as 




+ 


Ee,. 


, 1 , 

smmnA cOnt-\ — at 


Defining the following non-dimmsional parameters for the it th modal equation 


r* = Pk^ = non - dimensional time 


(Or, 


/q = — = imtial frequency ratio 
‘ Pk 


CO 


r = — = instantaneous frequency ratio 
Pk 


a 

= — 
Pk 


= nondimensional sweep rate 


= — — = dynamic magnification factor (or displacenent ratio) 


a. 


a 


■= static deflection 


equation (3.36) can be written as 


dT\ 


dH 1 


+'T Qn sin/MM,(rnr +— r„r^) 


Equation (3 38) being linear the individual respcaises for the force harmomcs 
calculated indqiendently and thai summed together to obtain the total response. 

The response for the zeroth harmonic is, readily 


(3 35) 


(3 36) 


(3.37) 


(3.38) 

can be 
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‘^OA'j 


(3.39) 


while for the higher harmonics, the response, as in the earlier chapter can be written using 
Duhamel’s Int^al, as (for zero initial conditions) 


= sin(J't -^■i)cosmn 

0 \ 2 

^(m+6)k (r) = Je sin(r* - r'*)sin»i«,[ 

0 V 


.,2 


dr' 


'ok 



dr' 

k 


(3.40) 


The nondimensional modal response from equations (3 39) and (3 40) above, can be 
processed along with equations (3 37), (3 33), (3 30) and (3 39) to obtam the blade displacement 
and moment fields Stresses are computed using Strength of Material formulae and the basic 
design stress is obtained by employmg the von-Mises critenon on pnndpal stresses 

3.3 Results 


A corr^uter program, based on the above formulations, is developed for obtainmg the 
forced vibration response of a turbme blade for step-up and stqi-down operations 

Fig 3.4 shows a typical turbine blade, used m the present study for illustration The 
geometnc and matenal properties of this blade are given in Table 3 1 . The first five natural 
frequenaes of the blade are listed m Table 3 2 Typical modal damping data and nozzle exatation 
for the blade have been shown in Figs (3.5) and (3.6) respectively These eiqienmentally generated 
data have been takoi from Vyas[1986] Fig 3 7 is the Can^bell diagram for the blade, which 
defines the resonant speeds from the interaction of the blade natural modes and the excitation force 
harmomcs (e g the interaction of the first blade mode with the first harmonic of the exatation 
force gives the resonant speed of 816 RPM, while the mterachon of the first blade mode with the 
sixth harmorac of the exatation force gives the resonant speed of 135 RPM). 

The program is illustrated for start 151 operation of the rotor from an initial speed of 0 
RPM to a final speed of 4000 RPM. The blade stresses have been computed for the resonances 
occurring, m this speed range, due to the mteraction between the first blade mode and first six 
harmomcs of the nozzle excitation. The computation is illustrated for typical accelerations of 2000 
RPM/min, 4000 RPM/min, 6000 RPM/min, 8000 RPM/min and 10000 RPM/min The effective 
nondimensional sweqi rates mrt^{alp^), for these accelerations, fall in the range 
1 X 10^ to 5 X 10'^ (with p^= Pi = 492 88 Hz , = 36, m = 1,2. . 6 ). 
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Fig. 3.4 


Blade profile. 
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TABLE 3.1. Geometric & Material Properties of the Blade 


Modulus of elastiaty {E) 

200 GPa 

Modulus of ngidity (G) 

72GPa 

Density ip) 

SOOOkg/m^ 

Blade length 

111 mm 

Disc radius 

180 mm 

Pretwist angle 

0 0 rads 

Stagger angle 

1 5708 rads 

Area of cross section (A ) 

17 8447x10'^ 

Moment of area about xr axis (/*,) 

38 39254x10'* 

Moment of area about yy axis (/_^.) 

80 08938x10'* /w^ 

Moment of area about xy axis {Jxy) 

46 33235x10'* /n^ 

X co-ordinate of centre of flexture 

: 27 45638 mm 

y co-ordinate of centre of flexture 

- 11.76078 mm 



TABLE 3,2. Natural Frequencies for first 5 Modes. 


MODE 

FREQUENCY (Hz) 

1 

492.9 

2 

1767.0 

3 

3067.9 

4 

4180.0 

5 

8586.5 








Fig. 3.7 
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RPM 


Campbell diagram. 
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TABLE 3.4 Transient Resonant Stresses 
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For such small effective swe^ rates, it is difficult to obtain a plot of the cyclic vanation 
of the stress for the entire sweep range, as solutions required to be computed for very large number 
of points on the range However, from the design point of view , the excessive stresses incurred by 
the blade while passing throu^ a resonant rotor speed would be of interest Therefore, for 
computational efficiency only the peak responses at the instances of momentary resonances are 
computed This is done, by employing the results obtained in nondimensional form for the case of 
a sin^e d^ee of freedom system Comparison of the modal equation (3.38) of the blade with 
equation (2 5) of the previous chapter and between the blade modal response equations (3 40) with 
equation (2 7) reveals, that the results for the single d^ee freedom system can be readily 
extrapolated to obtam the blade modal response Fig 3 8 shows a nondimensional instantaneous 
peak response surface, plotted from the results of Table 2 1 The nondimensional peak response is 
plotted, m this figure, as function of the two parameters, namely sweep rate and damping This 
surface has beai built into the computer program and the response peak response for the blade, for 
a given rate of acceleration, and gjvai dan^mg data, is directly read from the surface The blade 
stresses thus obtamed at momentary resonances durmg the start-up operation, are listed m Table 
3 3. It can be noted that the stress levels for lower sweep rates are higher, which is m consonance 
with the basic prmciples of the earlier chapter 
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3.4 Remarks 

The computer program has been illustrated only for a start-up operation However, it is 
capable of handling shut-down operations (i e non-zero initial speeds and n^ative rates of 
acceleration) The peak stress level expenenced by the blade, during a deceleration operation are 
not expected to be very different from the corresponding ones for positive acceleration, as can be 
seen from Table 2.1 of the previous diapter. The influence of initial speed on the peak response 
has not been, presently investigated However, the influence is not expected to be notable since the 
free vibration due to initial conditions is expected to wither away relatively fast. 



CHAPTER 4 


CONCLUSIONS 

An attarpt has beei made in the present study to estimate the resonant stresses 
momentarily experienced by a turbine blade during a variable speed operation like start-up. Ihe 
analytical formulation employs the Reissner’s functional in conjunction with Ritz process to wnte 
the equations of motion The eigenvalue problan is solved to obtain the free vibratiwi 
diaractenstics of the blade Modal analysis is used to process the equations of forced motion For 
a computationally effiaent calculation of the modal peak responses, the noRdimensional 
response data of the single degree freedom system is enqiloyed as a built-in data surfrce, 
in the blade program, whereby the blade modal responses are directly extrapolated from 
this surfrce The stress response of the blade is computed using Strength of Material formulae. 
The computer program developed on the basis of the analysis is illustrated for typical start-up 
operations. 

Calculation of such transient resonant stresses is relevant for a blade designed to work 
satisfactorily at an innocuous constant speed for power gaieration, may fail entirely due to the 
fatigue damage cumulated through rqneated infliction of high momentary resonant stresses, during 
every start-iqi or shut down operaration 


As a next stage of analysis, fatigue damage calculations can be taken up. Cumulative 
fatigue theories, like the Miner's Rule m conjunction with familiar S-N diagrams or more advanced 
fracture mechanics principles, taking the transient resonant stresses as inputs, can be err^iloyed to 
attempt blade life predictions. The present analysis employs a nondimensional peak response data 
surface for effiaent computation For fatigue life calculations, the envelope of the nondimensional 
response will be required to mark off time periods, during which the transient response is higher 
than some specified safe limit, say the endurance limit of the blade material. An approach to 
building such information withm a life estimation software would be to set up an empirical 
relationship for the response envelope by applying an appropriate cmrve fitting algorithm on the 
data generated in the presort study. 
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APPENDIX 


TTie equation of motion for the turbomadiine blade can be wntten in matnx form as 


^11 ^12 ^16 
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^ 2 \ ^22 ^26 
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: 1 

^ oosmut + 
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,Qo6^ 
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|Q ... . 

Os 

m 

Qm+6,6 


(A.1) 


where the sub matrices- 

I 

0 

/ 

Ms =p\^fifgdz 

0 

/ 

M22 =p^AfJgdz 

0 

I 

M^.,=p^Ar^gJjdz 

0 

/ 

=p\Af,rygj dz 
0 
/ 

=p\^rJ,gjCk 

0 

/ 

^33 = p\ [p^i^x +ry) + Jcg ]g, Sj dz 
0 



K^^=-pw^ 



+ sm^ 






-!{{''yS.+r,g;y,](m,+i^)dz 

0 



dz + ^sin2^j Ar^gJjdz 

0 0 


A^i = sm2 


<f>l4fzfjdz 
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K,,=-PW^: 

K23 = 
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-J f/fji^i + 4 )^ + cos^ <!)J AfJj dz 
■ 0 0 . 
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, 0 




/ 

K,,=\h,g'jdz 

0 

/ 

K,, = -\f;hjdz 

0 

K - \ 

^ ~ J W/2 - / / ^ 


dz 




e(i^ -/ / ^ 

'^V ^i>'i ^1 / 

J 

K,,=-\f:h^dz 


dz 


^54 - J 

0 

/ 


''.^^*1^1 


_ 7 7 \ 




"lW/2 -J J \ 

I 


dz 


K,,=\g:hjdz 


K 


rkh. 


66 


-F, 


-dz 


1 

I 

1 

Qo\ ^ ^^Qxfj ^ 

Qm.l =l^mJjdz 

2m+6,l ” J 

0 

o 

o 

1 

J 

1 

Qo2 ” J^0>»/; ^ 

Qm,2 

^m+6,2 ”” J* 

0 

0 

0 

/ 

1 

J 
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\b,^Sj 


dz 


All other elements are zero 


(A-2) 


{qf = {a„B,,C„D,] 


(A-3) 



